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The non-equilibrium structural and dynamical properties of semiflexible polymers 
confined to two dimensions are investigated by molecular dynamics simulations. 
Three different scenarios are considered: The force-extension relation of tethered 
polymers, the relaxation of an initially stretched semiflexible polymer, and semiflexi- 
ble polymers under shear flow. We find quantitative agreement with theoretical pre- 
dictions for the force-extension relation and the time dependence of the entropically 
contracting polymer. The semiflexible polymers under shear flow exhibit significant 
conformational changes at large shear rates, where less stiff polymers are extended 
by the flow, whereas rather stiff polymers are contracted. In addition, the polymers 
are aligned by the flow, thereby the two-dimensional semiflexible polymers behave 
similarly to flexible polymers in three dimensions. The tumbling times display a 
power-law dependence at high shear rate rates with an exponent comparable to the 
one of flexible polymers in three-dimensional systems. 
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I. INTRODUCTION 



Semiflexibility is a characteristic property of a broad range of biological polymers. Promi- 
nent examples are DNA, filamentous actin, microtubules, or viruses such as fd-virusesi~-. 
The rigidity is fundamental for their biological functions. For example, the DNA persistence 
length strongly affects its packing in the genome or inside a virus capsid. Actin filaments 
are an integral part of the cytoskeleton and their rigidity determines its particular mechan- 
ical properties. Hence, considerable effort has been devoted to unravel the mechanical and 
dynamical properties of semiflexible polymers^ - — . 

Advances in single-molecule spectroscopy prompted experimental and theoretical studies 
of non-equilibrium properties of semiflexible polymers^— . Fluorescence microscopy studies 
on single DNA molecules in shear flow reveal large conformational changes and an intriguing 
dynamics, denoted as tumbling motion^— . This implies specific non-equilibrium conforma- 
tional, dynamical, and rheological properties, which have been analyzed experimentally^— , 
theoretically^— , and by computer simulations^^—. 

These studies typically consider semiflexible polymers in three-dimensional space. Much 
less attention has been payed to polymers in two dimensions, although we may expect to 
see particular features in their equilibrium and non-equilibrium dynamical properties. Two- 
dimensional behavior is realized for strongly adsorbed polymers at, e.g., a solid surface, a 
membrane, or at the interface between immiscible fluids. Experiments reveal a strong de- 
pendence of the diffusive dynamics of adsorbed polymers on the underlaying substrate^ 1 ^. 
Moreover, theoretical and simulation studies predict a strongly correlated dynamics in two- 
dimensional polymer melts^. Little is known about the non-equilibrium properties of poly- 
mers in two dimensions. Here, we refer to recent simulation studies of end-tethered semi- 
flexible polymers, where the central monomer is periodically excited^ 1 ^. These simulations 
find a crossover from a limit cycle to an aperiodic dynamics with increasing stiffness. 

There are two major differences to three-dimensional systems. First of all, excluded 
volume interactions play a more pronounced role. The non-crossability leads, e.g., to a seg- 
regation of polymers in two dimensions^. We expect a strong impact of these interactions on 
non-equilibrium properties too. Secondly, hydrodynamic interactions can be neglected under 
certain circumstances^. This applies to strongly adsorbed polymers, where the polymer- 
substrate interaction dominates the dynamics of the polymer. It is certainly not appropriate 
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for polymers confined at fluid-fluid interfaces. 

In this article we investigate the non-equilibrium structural and dynamical properties 
of semiflexible polymers by computer simulations. As discussed above, we assume that 
the local polymer friction is determined by its interaction with the substrate and, hence, 
neglect hydrodynamics. Thus, we exploit the Brownian multiparticle collision (B-MPC) 



dynamics approach described in Refs. 



78-80 



By varying the chain stiffness, we gain insight 
into the dependence of the polymer properties on stiffness. Moreover, by comparison with 
existent results on three-dimensional systems, we uncover specific effects of the reduced 
dimensionality. 

Three different situations are considered. We briefly touch the force-extension relation 
of a semiflexible polymer and show that it is well described by theory*^. In addition, 
we examine the end-to-end vector relaxation behavior of initially stretched polymers. We 
find excellent agreement with the power-law dependence obtained in experiments^. This 
suggests that our model is a useful coarse-grained representation of a DNA molecule, at 
least for the considered properties. The major focus of the paper is on the non-equilibrium 
properties of semiflexible polymers under shear flow. We discuss a broad range of structural 
and dynamical quantities and stress the universal character and/or their particular, two- 
dimensional features. 

The paper is organized as follows. In Sec. [Ill the polymer model is described and the sim- 
ulation method is introduced. In Sec. II III the force-extension relation of a tethered polymer 
in a uniform external field is discussed. Section IIVI presents results on the time dependent 
relaxation behavior of a stretched semiflexible polymer. The structural properties of free 
polymers under shear flow are discussed in Sec. |V] and their tumbling dynamics is analyzed 
in Sec. I VI I Finally, Sec. IVHI summarizes our findings. 

II. MODEL AND METHOD 

The polymer is modeled as a linear chain composed of N beads of mass M. Its in- 
tramolecular interactions are described by the potential U = Ubond + Ub en d + U ex . Successive 
beads are linked by the harmonic bond potential 

JV-l 

U bond = y J2(\r i+1 - n\ - r ) 2 , (1) 



where is the position vector of bead i (i = 1, . . . , N), is the spring constant, and r 
the bond length. The bond bending potential 

JV-2 

Ubend = K^i 1 ~ COS (pi) (2) 
i=l 

accounts for the bending stiffness of the polymer, with k the bending rigidity and tpi the 
angle between two consecutive bond vectors. In the semiflexible limit k — >■ oo, the bending 
stiffness is related to the persistence length by L p = 2Kro/ksT, where k B T is the thermal 
energy, with T the temperature and k B Boltzmann's constant. Excluded- volume interactions 
are ensured by the shifted and truncated Lennard- Jones potential 



where r denotes the distance between two non-bonded beads and 0(r) is the Heaviside 
function (0(r) = for r < and 0(r) = 1 for r > 0). The dynamics of the beads 
is described by Newton's equations of motion, which we integrated by the velocity- Verlet 
algorithm with time step Atp^^. 

The polymer is coupled to a Brownian heat bath, which we implement via the B-MPC 
approach^ 1 ^^ 3 -. Hence, no hydrodynamic interactions are taken into account. In B-MPC, a 
bead performs stochastic collisions with a phantom particle which mimics a fluid element of 
a certain size. The momentum of the phantom particle is taken from a Maxwell-Boltzmann 
distribution of variance Mk B T and mean given by the average momentum of the fluid field, 
which is zero at rest or (Mjyi, 0) T in the case of an imposed shear flow of shear rate 7 in 
the xy— plane. For the stochastic process itself, we apply the stochastic rotation dynamics 
realization of the MPC method^^ 1 ^. Here, the relative velocity of a polymer bead, with 
respect to the center-of-mass velocity of the bead and the associated phantom particle, is 
rotated in the xy— plane by angles ±<x The time interval between collisional interactions is 
At, which is larger than the time step At p . 

The simulations are performed for the parameters a = 130°, At = 0.1t u , where the time 
unit is t u = ^mr%/(k B T), M = 5m, K h r 2 /(k B T) = 4 x 10 3 , e/(k B T) = 1, a = r , N = 51 
so that polymer length is L = 50ro, and At p = 10~ 2 At. With this choice for /%, the length 
of the polymer is kept constant within 1% for all systems. 




(3) 
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FIG. 1. Polymer extension along the direction of the external force for the persistence lengths 
Lp/L = 2(o), 4(«) as a function of the dimensionless force Fj = L p F(N — l)/(/c#T). The full line 
represents Eq. @ and has the slope —1/2. 



III. POLYMER FORCE-EXTENSION RELATION IN UNIFORM FIELD 



We consider a single tethered polymer with its endpoint T\ fixed at r\ = without 
additional restrictions on the orientation of the first bond. Every monomer is subjected 
to the external force F along the ^-direction of the Cartesian reference frame, e.g., due to 
an external electric field E x , where the force is F = qE x , with q the electric charge of the 
monomer. As is well known, the polymer is stretched along the force direction, where the 
extension increases non-linearly with increasing force^ 1 ^. Theoretical calculations based on 
the Kratky-Porocr^ wormlike chain model predict the asymptotic dependence 

; rji \ 1/2 

KrT 



L 



(4) 



2L P F(N- 1) 
for \xn\ —> L. 

Simulation results are presented in Fig. [T]for the persistence lengths L p /L = 2, 4. They 
agree very well with the theoretical prediction (J3J). This confirms that the model represents 
a continuous semiflexible polymer over the presented range of forces. Due to the discrete 
nature of the model, deviations will appear from the predictions of continuous semiflexible 



polymers 



for large forces, as shown in Refs. 



86 



90, and 



91 



, and a crossover will occur 



10 
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FIG. 2. Relaxation of the end-to-end distance along the stretching direction \xN(t/t u ) — xn(0)\ of a 
semiflexible polymer with L p /L = 2. The line indicates a fit in the time range t/t u = 10 3 — 50 x 10 3 
and has the slope 0.45. 



from the force-extension relation of a semiflexible model to that of a freely jointed chain. 
The effects of attractive interactions between non consecutive nearest neighbor beads^i on 
the described picture might be interesting to be investigated in the future. 



IV. RELAXATION OF STRETCHED POLYMER 

Releasing the force on a stretched polymer leads to its collapse into an equilibrium confor- 
mational state. This relaxation exhibits a characteristic time dependence. Figure [2] displays 
the relaxation behavior of the x-component of the end-to-end vector of a stretched semi- 
flexible polymer. The initial average stretching along the x— direction is xn{0)/L = 0.945, 
induced by the force FL p (N — 1)/ {ksT) = 200. We observe a power-law decrease of the ex- 
tension according to \xN(t) — xn(0)\/L ~ t 1 over a broad time scale. A fit of the data, which 
are averages over 20 independent realizations, yields the exponent 7 = 0.45 ± 0.03. This 
value is in remarkable agreement with the exponent 0.46 ± 0.08 found in two-dimensional 
experiments^. 
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V. SEMIFLEXIBLE POLYMER UNDER SHEAR FLOW 



In our study of semiflexible polymers under shear flow in two-dimensional space, we 
consider the persistence lengths L p /L = 0.1,0.4,2, 10. The corresponding equilibrium end- 
to-end vector relaxation times are r /t u ~ (161, 370, 676, and 707) xlO 3 , respectively. The 
strength of the flow is characterized by the Weissenberg number Wi = 7T0 in the range 
1< Wi < 800. 



A. Conformational Properties 

1. End-to-End Vector 

Probability distribution functions (PDFs) of the polymer end-to-end distance R e = \r^ — 
7*1 1 are presented in Fig. [3] for the Weissenberg numbers Wi = 8, 80, 800 and the various 
persistence lengths. As shown in Fig. [3](a), for small persistence lengths, the polymers are 
able to assume coil-like conformations, which give nearly constant PDFs over a wide range of 
end-to-end distances. Very small distances are suppressed by excluded-volume interactions 
and large distances are rarely sampled due to entropic penalties. However, larger shear 
rates lead to a sampling of large R e values. At low shear rates, an increasing persistence 
length naturally leads to a preference of large R e values. Shear, however, leads to an opposite 
behavior. At a given L v j L, an increasing shear rate gives rise to an increase in the probability 
distribution at smaller end-to-end distances. The effect becomes more pronounced for larger 
stiffnesses (cf. Fig. [3] (c) and (d)). This is in agreement with the predictions of Ref. 



49 



that 



semiflexible polymers under shear flow behave more and more like flexible polymers with 
increasing Weissenberg number. 

Figure H] displays the mean end-to-end distances as function of Wi. For every shear rate, 
(R e ) is smaller for the more flexible polymer. However, (R e ) increases for flexible polymers, 
whereas it decreases for the stiffer ones. As predicted by theory^, we expect that the end- 
to-end distances become similar for all stiffnesses in the asymptotic limit Wi — > 00. As the 
figure clearly reveals, in the stationary non-equilibrium state, a semiflexible polymer is never 
fully stretched. This has also been observed in simulations of flexible polymers^ 1 ^ and in 
experiments on DNA molecules^ 1 ^. 

Mean square end-to-end distances (R 2 ex ) along the flow direction are displayed in Fig. O 
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FIG. 3. Probability distributions of the polymer end-to-end distance R e = \r^ — n| for the 
Weissenberg numbers Wi = 8(o), 80(«), and 800(*) and L p /L = 0.1(a), 0.4(6), 2(c), 10(d). 



They also show that semiflexible polymers are never fully stretched. Moreover, the various 
curves reveal a very weak persistence length dependence for polymers with L p /L > 0.4. 
They closely follow the same Weissenberg number dependence. This has been predicted in 
Ref. |49| and is related to the fact that the rather stiff polymers become first aligned with 
the flow at moderate shear rates. Only at larger shear rates, deformation sets in. This is 
also evident from Fig. HJ which clearly exhibits a dependence of the "critical" Weissenberg 
number on the persistence length, above which (R e ) decreases with increasing shear rate. 
Below the critical value, the polymers are aligned by the flow and above, in addition, they 
are deformed^. 

In Fig. E]mean square end-to-end distances {R^ y ) are shown along the gradient direction. 
The polymers of the various stiffnesses shrink transverse to the flow direction. Thereby, 
we observe a slight dependence on persistence length over the considered range. The decay 
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FIG. 4. Mean values (R e ) of the end-to-end distance as function of the Weissenberg number for 
the persistence lengths L p /L = 0.1(o), 0.4(»), 2(*), and 10(A). 



at larger Wi can approximately be described by the power-law {R 2 ey ) ~ Wi~ u , with v fa 
1/2. This dependence is consistent with the decay of the radius-of- gyration tensor of three- 
dimensional systems2^2^. However, theoretically an exponent v = 2/3 is^expectecr^ 1 ^ 
which seems to be reached for much higher Weissenberg numbers in Ref. 
exponent v — 1/2 could characterize a crossover behavior only. 
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Hence, the 



2. Bond Angle 

To further characterize the polymer conformational properties, we present in Fig. [7] the 
average bond angles (tpi) (f2]) between successive bond vectors along the semiflexible poly- 
mers. As expected, the angles (ifi) decrease with increasing persistence length and are close 
to zero for L p /L = 10. At small persistence lengths, the values of (ipi) decrease with increas- 
ing shear rate, specifically toward the middle of the chain, due to polymer stretching by the 
flow. The situation is reverted at large persistence lengths, where the angles (cpi) increase 
with the shear rate. 
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FIG. 5. Mean square end-to-end distances along the flow direction as function of the Weissenberg 
number for the persistence lengths L p /L = 0.1(o), 0.4(»), 2(*), and 10(A). \R^ ) is the mean 
square end-to-end distance at equilibrium. 



Average bond angles 

1 N-2 

^ = J^H fo> ( 5 ) 

i=i 

are displayed in Fig.[8]as function of Weissenberg number and for various persistence lengths. 
Evidently, the mean values are independent of shear rate for the larger persistence lengths. 
Only for the considered most flexible polymer a decrease of ((p) is found as already expected 
from Fig. [7J 

This minor change in the bond-bond orientational behavior is surprising in the light of 
the decreasing mean end-to-end distance (cf. Fig. HJ). This is explained on the one hand by 
the nearly rigid rod-like rotation of the semiflexible polymers at lower Weissenberg numbers 
and on the other hand by the formation of U-shaped conformations with only small and local 
bending of the polymer, as reflected in Fig. [7Jc) and (d) at higher values of Wi. Typical 
conformations of stiff polymers at low and high Weissenberg numbers are presented in Fig. [9] 
(see also the movies in the supplementary material^). 
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FIG. 6. Mean square end-to-end distances along the gradient direction as function of the Weis- 
senberg number for the persistence lengths L p /L = 0.1(o), 0.4(»), 2(+), and 10(A). The full line 
has the slope —1/2. (R^o) is the mean square end-to-end distance at equilibrium. 



B. Alignment 

Polymers under flow are not only deformed, but also exhibit a preferred, flow induced 
orientation^'^'^'^^i. To characterize the degree of alignment, we calculate the prob- 
ability distribution of the angle between the end-to-end vector and the flow direction. 
Examples are shown in Fig. [10] for various Weissenberg numbers and persistence lengths. 
There is no preferred angle at equilibrium. With increasing shear rate, the distribution func- 
tion exhibits a maximum at a non-zero, positive value <p m . This maximum shifts to smaller 
values with increasing shear rate. At the same time the distribution function becomes nar- 
rower. The latter implies that a polymer aligns preferentially in a particular direction and 
samples other angles only rarely. 

Interestingly, the probability distribution functions become asymmetric with increasing 
shear rate and exhibit a second maximum at large angles. Thereby, the second maxi- 
mum moves to larger angles with increasing stiffness. The strong asymmetry, particular 
for lower Weissenberg numbers, seems to be specific for two-dimensional systems, because 
(flexible) polymers in three dimensions^ exhibit more symmetric distributions. The second 
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FIG. 7. Average local bond angle {(fi) along the polymer contour for the Weissenberg numbers 
Wi = 8(o), 80, (•), and 800(*) and the persistence lengths L p /L = 0.1(a), 0.4(6), 2(c), and 10(d). 



peak is peculiar for rather stiff polymers. However, it is not clear whether it appears in 
two-dimensional systems only. Here, studies of three-dimensional semiflexible polymers are 
required to resolve the issue. We report for completeness the fact that the appearance of 
two peaks in P{4>) at finite values, symmetric with respect to <p — 0, was observed for two- 
dimensional grafted polymers with the first bond fixed along the x-direction and L p ~ L 
under equilibrium conditions.— 

As shown in Fig. (TTJ the probability distribution functions depend only weakly on the 
persistence length for Weissenberg numbers Wi < 10 2 . In particular, the position of the 
maximum is virtually independent of L p . 

Figure [T^l displays the angles <p m of the central maximum of P{4>)- For small Weissenberg 
numbers Wi < 10, tan(20 m ) decreases as Wi' 1 , whereas for larger Weissenberg numbers 
the dependence tan(20 m ) ~ Wi' 1 / 3 is obtained. A similar dependence is found for flexible 
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FIG. 8. Average bond angles {ip) ([5]) as function of the Weissenberg number for the persistence 
lengths L p /L = 0.1(o), 0.4(»), 2(*), and 10(A). 



polymers in three dimensions^ and is predicted theoretically^ 1 ^ independent of dimension. 
As already suggested by Fig. [TT], nearly the same degree of alignment is obtained independent 
of stiffnesses. However, the values of tan(20 m ) are somewhat smaller for the more flexible 
polymers in the range Wi > 10, as predicted theoretically^. 

Similarly, the width A0, which is defined as the full width at half maximum of the 
distribution function P((j>), decreases as A0 ~ Wi' 1 / 3 at high shear rates. This has also 
been observed in experiments^ and predicted theoretically^. 



VI. TUMBLING DYNAMICS 

As mentioned in the Introduction, polymers in shear flow undergo a tumbling motion. A 
characteristic tumbling time can be obtained from the distribution function P(t) of times 
between successive zeros of the end-to-end vector component R ex {t) along the flow direction. 
This distribution exhibits the exponential decay P(t) ~ exp (— t/r^) at large times, from 
which the tumbling time is extracted. Normalized tumbling frequencies ~ l/r</, are 
depicted in Fig. [131 where the full line has slope 2/3. 

As for a three-dimensional system, we obtain ~ Wi~ 2 ^ 3 for the shear rate dependence 
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FIG. 9. Typical conformations at consecutive times (from top to bottom) of a semiflexible polymer 
{L p /L = 10) at Wi = 8 (left panel) and Wi = 800 (right panel). 

of the tumbling times. This confirms that the tumbling times of semiflexible polymers exhibit 
the same asymptotic Weissenberg number dependence as flexible polymers. Moreover, two- 
dimensional tumbling behavior seems to be similar to three-dimensional one for semiflexible 
polymers. 

Interestingly, the frequencies for the larger persistence lengths seem to exhibit a fast, 
rather abrupt increase in the vicinity of Wi ~ 10 and approach the dependence Wi 2 ^ 3 
for large Weissenberg numbers only. Such a dependence has not been observed in three 
dimensions so far nor be predicted theoretically. Whether this is a specific two-dimensional 
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FIG. 10. Probability distributions of the angle (ft for the Weissenberg numbers Wi = 8(o), 80(»), 
and 800(*) and the persistence lengths L p /L = 0.1(a), 0.4(6), 2(c), and 10(d). 
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FIG. 12. (a) Angles (j) m of the maximum and (b) width A(p of the distribution function P((f>) as 
function of the Weissenberg number for the persistence lengths L p /L = 0.1(o), 0.4(»), 2(*), and 
10(A). The slopes of the full lines are —1 and —1/3, respectively. 



features needs to be addressed by simulations of semiflexible polymers in three dimensions. 
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FIG. 13. Normalized tumbling frequencies to/t^ as function of the Weissenberg number for the 
persistence lengths L p /L = 0.1(o), 0.4(»), 2(*), and 10(A). The line has the slope 2/3. 

VII. CONCLUSIONS 



We have presented results for the non-equilibrium structural and dynamical properties 
of semiflexible polymers confined to two dimensions. The analysis of the force-extension 
relation of a semiflexible polymer in a uniform external field yields excellent agreement 
with theoretical predictions^ 1 ^ and confirms that the applied model is very well suited to 
describe semiflexible polymers. Our studies of the end-to-end distance relaxation behavior 
of initially stretched semiflexible polymers confirm the experimentally obtained power-law 
time dependence t 1 with the exponent 7 = 0.45, which is in close agreement with the scaling 



prediction 7 = 1/2 of Ref. 



77 



We have also studied the conformation properties of semiflexible polymers under shear 
flow. We clearly find strong shear-induced conformational changes. Beyond a stiffness de- 
pendent Weissenberg number, the average polymer extension decreases with increasing shear 
rate, in contrast to flexible polymers, where the extension increases. Visual inspection shows 
that U-shaped conformations appear. Such conformations have also been observed for semi- 
flexible polymers in microchannel flows, both experimentally^ and in simulations^"^. This 
confirms that the end-to-end distances become similar in the asymptotic limit of infinite 
shear rate independent of the stiffness^. As for flexible polymers in three dimensions, the 
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semiflexible polymers preferentially align along the flow direction. However, the distribution 
functions of the end-to-end vector alignment angle are clearly more asymmetric at low Weis- 
senberg numbers than for flexible polymers and exhibit a second peak for large stiffnesses. It 
is not evident whether these effects are caused by stiffness or confinement to two dimensions. 
Further simulation studies are necessary to resolve this question. 

We have also shown that semiflexible polymers exhibit a tumbling motion, where 
the tumbling times approximately show the dependence ~ Wi~ 2 ^ 3 on shear rate. 
Hence, semiflexible polymers reveal in essence the same tumbling behavior as flexible 
polymera22^LasmZL21 and rods^a. 

There are various aspects, e.g., the appearance of a second peak in the orientational 
distribution functions and their asymmetry at low shear rates, which need further investi- 
gations to clarify the underlaying mechanism. This requires theoretical calculations and/or 
simulations in two and three dimensions. We hope that our results will stimulate such the- 
oretical studies as well as experimental investigations and will be valuable in the respective 
endeavors. 
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